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ABSTRACT 

Large-scale radial transport of solids appears to be a fundamental consequence 
of protoplanetary disk evolution based on the presence of high temperature min- 
erals in comets and the outer regions of protoplanetary disks around other stars. 
Further, inward transport of solids from the outer regions of the solar nebula 
has been postulated to be the manner in which short-lived radionuclides were 
introduced to the terrestrial planet region and the cause of the variations in oxy- 
gen isotope ratios seen in primitive materials. Here, both outward and inward 
transport of solids are investigated in the context of a two-dimensional, viscously 
evolving protoplanetary disk. The dynamics of solids are investigated to deter- 
mine how they depend on particle size and the particular stage of protoplanetary 
disk evolution, corresponding to different rates of mass transport. It is found 
that the outward flows that arise around the disk midplane of a protoplanetary 
disk aid in the outward transport of solids up to the size of CAIs and can increase 
the crystallinity fraction of sihcate dust at 10 AU around a solar mass star to as 
much as ~40% in the case of rapidly evolving disks, decreasing as the accretion 
rate onto the star slows. High velocity, inward flows along the disk surface aid 
in the rapid transport of solids from the outer disk to the inner disk, particu- 
larly for small dust. Despite the diffusion that occurs throughout the disk, the 
large-scale, meridonal flows associated with mass transport prevent complete ho- 
mogenization of the disk, allowing compositional gradients to develop that vary 
in intensity for a timescale of one million years. The variations in the rates and 
the preferred direction of radial transport with height above the disk midplane 
thus have important implications for the dynamics and chemical evolution of 
primitive materials. 
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Introduction 



The results of the Stardust mission, which collected roughly 0.001 g of material ejected 



from Comet Wild 2 (Brownlee et al. 2006), suggests that large-scale transport in the solar 



nebula was a fundamental process in shaping the properties of the planetesimals that formed 
within it, and by extension, the planets that accreted from these building blocks. Among 
the materials returned to Earth were crystalline Mg-rich silicate grains which require tem- 
peratures in excess of 1000 K to form from the amorphous precursors expected to have been 



inherited from the natal molecular cloud (Hallenbeck et al. 2000), or even higher tempera- 



tures to condense directly from a gas phase (Wooden et al. 2007, and references therein). 
Further, a refractory grain called Inti, whose mineralogy and oxygen isotopic compositions 
is reminiscent of the Calcium, Aluminum-rich Inclusions (CAIs) found within chondritic me- 



teorites was also found within the Stardust aerogel (Brownlee et al. 2006 Zolensky et al 



2006 McKeegan et al. 2006). These findings have led many to conclude that large amounts 



of materials that formed or were processed in the hot, inner regions of the solar nebula mi- 
grated outwards tens of AU to be accreted by comets. The exact distance that materials 
were transported is still the matter of debate as comets likely originate from beyond 20 AU 
( Charnoz and Morbidelli||2007 ), but the cometesimals from which these objects formed likely 
originated anywhere from ~5 AU (Supulver and Lin 2000) to ~100 AU (Weidenschilling 



1997, 2004). 



That dust was mobile in the solar nebula is not surprising, as we have long recognized 
that transport processes must have operated in our solar nebula. For example, the chondritic 
meteorites, which are relatively unaltered products of our solar nebula, come in a variety of 
types, each characterized by their bulk chemistry, oxygen isotopic ratios, as well as chondrule 
sizes and abundances (see review by Weisberg et a/.|[2006 ). That each of these bodies record 
such different formation environments despite forming in close proximity to one another (a 
scale of 1-3 AU in the solar nebula, in all likelihood) and over a time period of 1-3 million 
years (e.g. Kita et al]|2005 ) suggests that the solar nebula was either spatially heterogeneous 
or that its chemical and isotopic properties changed with time due to the transport of 



materials within (Stevenson and Lunine 1988 Cyr et al. 1998 Cuzzi et al. 2005; Ciesla and 



Cuzzi 2006) 



Further, the components of chondritic meteorites-the refractory inclusions, chondrules, 
and matrix which contains pre-solar grains-all record very different chemical and isotopic 
environments as well as different thermal histories prior to their incorporation into common 
meteorite parent bodies. It is difficult to imagine how such diverse materials could be brought 
together without invoking solar nebular transport of some kind. 



Fortunately, we can gain insight into the details and driving mechanisms for transport 
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in our solar nebula by examining other protoplanetary disks. Crystalline silicates are abun- 
dant in a number of disks around other stars at levels greater than those seen in molecular 
clouds, particularly in the hot inner regions, indicating that their formation occurs at high 
temperatures, while their presence in the cool outer regions of disks is taken as evidence 



that large-scale radial transport of these materials occurred 
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). That evidence for this 



processing and transport is seen in disks around stars of all masses suggests that this is a fun- 
damental consequence of disk evolution. Thus by understanding the processes that control 
or operate during protoplanetary disk evolution we may identify what processes controlled 
how solids were transported in our own solar nebula. 

Protoplanetary disks are known to be dynamic objects that transport mass inward to 
be accreted by their central stars as part of their stars' last stages of pre-main sequence 
evolution (see review by DuUemond et a/.||2007 , and references therein). Small dust particles 
that are entrained in the gas would also be subjected to these motions, thus providing a means 
of large-scale movement. The exact driving force of this evolution remains the subject of 
ongoing work, but likely involves instabilities of some kind, such as the magnetorotational 



instability (MRI) of Balbus and Hawley (1991) or gravitational torques arising from massive 
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Hartmann et al. (1998) demonstrated that the observed accretion rates of T-Tauri stars 



can be matched by adopting the classical a-viscosity model in which the evolution of the 



disk is treated as arising due to viscous stresses in the differentially rotating gas (Shakura 
and Sunyaev||1973 Lynden-Bell and Pringle||1974 ) . These stresses serve to drive mass inward 
with time and allow the disk to spread in the radial direction to conserve angular momentum. 
The a-viscosity model is a turbulent viscosity model as it assumes that turbulence within the 
disk acts as the source of the viscosity, because the molecular viscosity of the gas is orders 
of magnitude too small to explain the observed mass accretion rates. A similar model, the 
/^-viscosity model, has also been proposed in which the turbulence specifically arises due to 



hydrodynamic instabilities in the disk (Richard and Zahn 1999). This latter model assigns 



a viscosity to the disk whose functional dependence on disk properties differs slightly from 
the a-viscosity model. 

Because of the turbulent nature of the gas, the random motions associated with its evo- 
lution lead to a diffusivity of material within it, allowing dust particles to be redistributed in 
a manner to smooth out concentration gradients. This diffusion has previously been investi- 
gated to evaluate whether particles that formed or were processed in the hot, inner regions of 
the solar nebula could have been carried outwards against the inward flows associated with 



mass transport to be delivered to where comets formed ( Gail 2001 , 2004 Bockelee-Morvan 



et al. 2002 ) . In the case of Bockelee-Morvan et al. ( 2002 ) , they found that such a scenario 



was possible, however, in order to get a large fraction of the silicate grains in the outer nebula 
to be crystalline, in agreement with the observations of comets Halley and Hale-Bopp, the 
nebula must have started out "warm" and compact. That is, the entire nebula mass was 
initially distributed within 10 AU of the Sun, resulting in mass accretion rates in excess of 
10~^ M0/year. Such a mass accretion rate is greater than observed for disks around roughly 



solar mass stars (e.g. Calvet et al. 2005), though the data is lacking for ages less than 10^ 



years old. Gail (2001, 2004) also found that high mass accretion rates favored the outward 



transport of high temperature materials, though the fractions of crystalline forsterite in these 
models were small (~1%). 

Previous studies of dust transport in protoplanetary disks, such as those described above, 
have largely relied on one- dimensional models where the changes in the surface (or column) 
density of the dust are calculated as a function of time and location as the disks evolved. 
Thus dust at a given location of a disk is expected to follow the same dynamical behavior 
wherever it is located in that column. However, as discussed in detail in the next section, 
the dynamics of dust particles depend on the local gas density, pressure, and their respective 
gradients, all of which vary with height in a protoplanetary disk. The importance of such 
considerations for the chemical evolution of protoplanetary disks has been recognized in 



recent work by Tscharnuter and Gail (2007), Ciesla (2007b), and Wehrstedt and Gail (2008) 



in which the rates as well as the preferred direction of radial transport were found to vary 



with height above the disk midplane. Here, the model of Ciesla (2007b) is extended to 



examine the roles that these variations played in the transport of primitive materials in our 



solar nebula as well as other protoplanetary disks. The model of Ciesla (2007b) is used as it 



accounts for the motions that develop as a result of gas drag and vertical settling, which are 



not included in Tscharnuter and Gail (2007) and Wehrstedt and Gail (2008), and therefore 



can be applied to cases beyond vapor species and sub-micron dust particles. 

In the next section, the two-dimensional (radial and vertical) model of particle trans- 
port in protoplanetary disks is described, with the key differences with one-dimensional 
(radial) models highlighted. Section 3 describes the protoplanetary disks to which this two- 
dimensional model is applied, with the specific numeric approaches described in Section 4. 
Section 5 presents the results for transport calculations for disks of the various structures 
described. The implications that these results have for the evolution of primitive materials 
in our solar nebula and other protoplanetary disks are then discussed. 
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2. Model Description 



A key parameter in quantifying how a solid particle responds to the motions of the gas 
in a protoplanetary disk is given by the particle stopping time: 



where p is the material density of the solid particle, c is the local speed of sound, pg is the 
local gas density and a is the particle radius. This is approximately how long it takes for a 
particle to lose its relative velocity with respect to the gas in which it is embedded. Particles 
with short stopping times are very well coupled to the gas, moving almost as gas molecules 
themselves, whereas particles with long stopping times decouple from the gas and begin to 
follow a separate evolutionary path. Strictly speaking, this particular formulation is only 
true for particles whose radii are less than the mean free path of the gas. As the focus here 
is on the primitive materials found in chondritic meteorites, comets, or observed in disks 
around other stars, all of which arc no more than 1 cm in size, this particular formula can 
safely be used for all cases considered here. 

Given that the stopping time is inversely proportional to the gas density, it is not 
surprising that the particular dynamics of the particles will vary with location within a 
protoplanetary disk. Particles will generally have shorter stopping times and be trapped 
to the gas in the inner disk, whereas particles in the outer disk respond more slowly, and 
decouple to greater degrees. Similar variations in behavior are expected as one considers the 
gas density profile above the disk midplane which is given by: 



where pg^ is the density of the gas at the disk midplane, z is the height above the midplane, 
and H is the local gas scale height, given by c/VLk where VLk is the local Keplerian frequency. 
As the stopping time simply scales as the inverse of the gas density, this means that the 
variations with height will go as: 



In other words, the stopping time of a particle at z—H will be 2.7 times what it is at the 
disk midplane, at z—2H, it will be 55 times, and at z—3H it will be ~8000 times. Thus 
particles will respond differently to gas motions depending on where they are in the column 
of gas at a given radial location in a protoplanetary disk. 



Pg (z) = Pg,oexp {-Z'^/H^) 



(2) 




(3) 



Not only does the response of a particle vary with height, but the particular forces that 
it will be subjected to will vary with height as well. For example, gas drag migration of 
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particles arises due to the differential rotation of the gas and solids in a protoplanetary disk. 
This differential rotation is due to the fact that the solids feel a central force due to the 
gravity of the central star whereas the central force that the gas experiences is modified 



by the radial pressure gradient (Adachi et al. 1976 Weidenschilling 1977). Around the 



midplane, the radial pressure gradient is largely negative as disks are expected to have hot, 
dense gas near the star and cool, sparse gas further out. However, at higher altitudes, the 
radial pressure gradient will actually become positive. This is due to the fact that while the 
surface density of the disk generally decreases with distance from the star, its vertical extent 
grows with distance, except under the most extreme conditions (that is since H = c/Qk, as 
long as c does not fall off more steeply than r~2, or T as r~^, the scale-height will increase 



with distance from the star). Takeuchi and Lin (2002) demonstrated that for a vertically 



isothermal disk, whose surface density and temperature can be described by the power laws: 



S(r) 



/ r \P 



T(r) 



(4) 
(5) 



iw)' 

where r is the radial coordinate, and the AU subscripts denote the values of the surface 
density and temperature at 1 AU, the radial pressure gradient could be expressed by: 
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(6) 



With the pressure gradient being negative at the midplane, it switches sign at z= ± [(2p + q — S) / {q 
meaning that above this height the gas will have an orbital velocity which exceeds the Kep- 
lerian rotation rate. This produces a tailwind which pushes materials outwards for the same 
reasons that materials migrate inwards due to gas drag around the disk midplane. Thus the 
rate and preferred direction of transport via gas drag will vary with height above the disk 
midplane. 

The variations in transport become even more pronounced when one considers the mo- 
tions of materials in a viscously evolving disk. In one-dimensional (ID) disk models of viscous 
protoplanetary disks, the evolution of the gas is found by vertically integrating the stresses 
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exerted on each annuli of gas in the disk and determining the net radial motions (Shakura 



and Sunyaev 1973 Lynden-Bell and Pringle 1974). In reality, however, the stresses exerted 



within the disk will be functions of the local gas properties and will vary with height as these 
properties change, potentially leading to very different behavior than the net radial motions. 



These variations were investigated in previous studies by Urpin (1984), Takeuchi and Lin 



(2002), and Keller and Gail (2004) who found analytic expressions for the variations in the 



gas motions in the disk with height above the disk midplane. These authors found that, in 
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agreement with numerical simulations ( Kley and Lin|[l992 Rozyczka et a/.|[l994 ), meridonal 
flows develop around the disk midplane. That is, in steady-state disks in which the mass 
accretion rate of material in the disk is constant in time and location in the disk, the net-flux 
of gas continues to be inward. However, the inward movement of material is largely relegated 
to the gas in the upper part of the disk, whereas the gas around the midplane moved out- 
wards. This outward movement of material is due to angular momentum conservation-that 
is, as mass is moving inward through a disk, some amount is pushed outwards to conserve 
angular momentum. That this happens around the disk midplane is due to how the mass of 
the disk is distributed, in that the gas density has its steepest radial gradient along the disk 
midplane, resulting in stresses that drive material outward. 



The structure of the disk flow can be found, following Urpin (1984) .Takeuchi and Lin 



(2002) and Keller and Gail (2004), by solving the azimuthal component of the momentum 



equation for rotating viscous fluid in cylindrical coordinates: 
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simultaneously with the equation for mass conservation: 

d 



dp„ Id, , 
-m^rd-r^'P^'''^ 







where Vr and Vz are the radial and vertical velocities of the gas, Qg is the orbital frequency 
of the gas, and u is the viscosity of the gas. 

Thus, given a steady-state protoplanetary disk with a given structure, it is possible to 
calculate the rates of transport of particles at any location due to the large-scale viscous 
flows, gas drag migration, and turbulent motions. Calculating the dynamical evolution of 
particles in the disk requires solving the two-dimensional advection-diffusion equation in the 
radial and vertical directions: 



— + ——(rV C) + ——(V C)-— — 
dt Pgr dr ' ^ pgdz ^ pgr dr 



dC\ 1 d 



dC 

dz 



(9) 



where is the total radial velocity of the dust particles at a given location due to both the 
large-scale flows of the disk and gas-drag migration, Vz is the total vertical velocity of the 
particles due to both the large-scale flows and vertical settling due to gravity, and C is the 
concentration of the particles deflned as pa/ Pg where pa is the mass density of the particles 
in the disk. Thus the concentration, C, purely reflects the concentration of material relative 
to the ambient gas and not necessarily relative to other dust that would be present in the 
disk. 
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For this study, rather than assuming that the protoplanetary disk mass density and 
thermal structure can be defined by the power laws in Equations (4,5), a self-consistent 
solution for the disk structure is found given a particular mass accretion rate. The process 
by which the surface density and temperature are found at each location is described below. 
The vertical structure of the disk is found by assuming hydrostatic equilibrium. This is 
a bit of a simplification as the vertical motions of the gas imply that perfect hydrostatic 
equilibrium is not completely achieved. However, because of the relatively small values of 
the gas vertical velocities, both Takeuchi and Lin (2002) and Keller and Gail (2004) showed 
that this assumption is accurate to order if^/r^ (0.1 to 1%). 



3. Disk Structure 



The structure of the disks in the models considered here were found by assuming that 
the disk was in steady state, that is, that the accretion rate through the disk M, is constant 
with time and location. Under such an assumption, the surface density of the disk and the 
viscosity at a given location. 



r, are related by (|Lynden-Bell and Pringle 1974|): 
1 



Ez/ 



-M 1 




(10) 



where E is the local surface density of the disk and rj„ is the inner edge of the disk. 



The midplane temperature of the disk, T^, which helps define the viscosity, is determined 
by balancing the rates of heating with the rate at which the disk cools from radiation. Two 
sources of heating were considered here: that due to viscous dissipation and that due to 



irradiation from the central star. The rate of viscous heating in the disk is given by (Ruden 



and Pollack 


1991 


Stepinski 


1998) 
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K 



which under steady-state assumptions becomes: 

3 



47r 



1 - 




The heating rate due to irradiation from the star is given by: 

O = - ^* 



:iii 



(12) 



(13) 



where (p is the angle at which the photons from the star impact the disk surface and L^, is 



the luminosity of the star at the time of interest (Chiang and Goldreich 1997; Brauer et al. 
2008|. 
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Given that the coohng rate of the disk is given by: 

Qcooi = 2crTg^ 



(14) 



where Te is the effective temperature of the disk surface, the effective temperature of the 
disk due to each individual effect can be calculated as: 



SGMqM 
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respectively. These expressions give the effective temperatures of the disk if viscous dissipa- 
tion or stellar irradiation dominated the heating of the disk respectively. Here, however, the 
important quantity for determining the viscosity of the disk is the midplane temperature. 
This can be estimated using the expression: 

3 



) e,visc ' e,trr 



(17) 



where r is the optical depth from the disk surface to the disk midplane (C. Dullemond, pers. 
communication). This formula accounts for the the fact that the energy generated interior 
to the disk by viscous dissipation must diffuse to the surface before being radiated away. 
The optical depth is given by /tS/2, where k is the opacity of the disk. 

Given values of M, the stellar luminosity, and values for k, and u everywhere, the surface 
density of the disk can thus be solved using an iterative approach. Here, values of M=10~^, 
10~^, and 10~^ Mq/jt were investigated and are discussed below. In all cases, the stellar 



luminosity was given by L^=4:nR^aT^, with i?*=3i?0 and T*=4000 K as in Ruden and 



Pollack (1991) and Ciesla and Cuzzi (2006). The incidence angle was taken as 0=0.05 for 



all simulations here, a value typical for disks where H/r increases with distance from the 
central star (e.g. Brauer et a/.||2008 ). While the actual value of will vary with location and 
time in a disk, given that the heating rate goes as 04, slight changes in its value will not 
produce significant shifts in the final temperature. A simple opacity law was adopted, where 
K= 5 cm^/g everywhere ( Cassen|[l994 ). Figure 1 shows the surface density and temperature 
profiles for the disks explored here. 



4. Numerical Approach 



In order to solve for the dynamical evolution of particles in the disk, the structure of 
the disk is defined on a two-dimensional, cylindrical (r — z) grid. The grid is logarithmically 
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spaced in both the radial and vertical directions with rj_|_i/rj=1.05 and Zj^i/ Zj=l.l, where 
ro=0.5 AU and ^0=8^10^'' cm. Given the surface density and thermal structures of the disk, 
the density at each location on the grid is defined based on the assumptions described in the 
previous section. 

The large-scale flows of the gas are calculated using an iterative finite-difference scheme. 
The radial velocity of the disk is first calculated using equation (7) and assuming f ^=0. These 
values of the radial velocity are then used in equation (8) to determine the values of Vz- These 
updated values are then applied to equation (7) and then the process repeated until the 
velocities converge to values that do not change by more than 1% in subsequent iterations. 
An example of the calculated flow structure is shown in Figure 2, which is qualitatively 
representative of the flow fields throughout the disk for each mass accretion rate shown. 

As the particles to be considered here are largely small with ts much less than an orbital 
period throughout the protoplanetary disk, their motions due to the large-scale motions of 
the gas are assumed to be the same as that of the gas. The velocities of the particles due 



to vertical settling and gas drag migration are calculated as in Tanaka et al. (2005). The 
net velocities of the dust particles at any point in the disk are thus found by summing those 
due to the large-scale flows of the disk and those due to gravitational settling and gas drag. 
Figure 3 shows the outward flow rates of the gas and the inward migration rates of the a=0.5 
mm particles at the disk midplane for the different disks considered here. For comparison. 



the inward drift rates of the a=5 /im dust would be ~100x smaller at the midplane (Cuzzi 



and Weidenschilling 2006). 



All transport calculations are done using a van Leer (1977) interpolation scheme. This 
significantly limits the amount of numerical (or artificial) diffusion that would occur as the 
calculations are performed. The calculations are explicit in nature, and thus the timestep is 
limited to the smallest time it takes for materials to advect or diffuse across a grid cell, thus 



satisfying the Courant-Friedrichs-Lewis stability criterian (Press et al. 1992). 



In order to ensure stability in the calculations, all transport in the disk is limited to 
occur within three gas pressure scale-heights of the disk midplane. It is within this range 
that more than 99% of the mass of the disk is contained, allowing such a consideration to 
be made without significant loss of accuracy. To impose this condition, all vertical velocities 
above z=3H at all radial locations are set to be equal to -10 km/s, and the diffusivity of all 
particles are set to zero. This forces all particles that are transported to these high altitudes, 
either through vertical diffusion or through radial transport from an adjacent radial bin 
which is below the 3H computational roof, to rain back down to the computational grid 
without any loss of mass from the simulation. 
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This approach contrasts with the two-dimensional transport models of [Keller and GaiT 



(2004) and Wehrstedt and Gail (2008) which used spherical geometries to study disks that 



extended 10° and 4.01° above the disk midplane respecitvely. While a direct translation into 
the effective scale-heights that these figures correspond to is difficult as it would depend on 
the temperature and location in the nebula, they would be roughly ~l-2if for typical nebular 
parameters. For comparison, the scale-heights at 1 and 20 AU for the different disk models 
considered here are listed in Table 2. These heights correspond to 3-8° above the midplane 



at 1 AU and 4-8° at 20 AU. Thus Keller and Gail (2004) and Wehrstedt and Gail (2008) 



focused on the region of the disk below roughly one scale height, and did not consider the 



motions that developed above this height (see, for example. Figures 3 and 4 of Wehrstedt and 



Gail||2008 ). Given the flow structure of the disk, this means that the rapid inward transport 
of the gas at high altitudes (above ~ 0.35-0.5 AU in Figure 2) are not modeled. While the 
gas density (and under well-mixed conditions the solid density as well) at these altitudes is 
much less than around the midplane, the rapid inward velocities at these altitudes means 
that the mass flux is still significant. Neglecting these higher altitudes leads to focusing on 
the region of the disk where outward transport of materials is more likely, and thus would 
lead to overestimates of the outward fluxes of materials in a protoplanetary disk. 

The computational domain was chosen to ensure that its vertical scale exceeded the 
computational roof throughout the disk in the model. While this roof makes the particular 
choice of boundary conditions on the upper level unimportant, dC jdz = was used in the 
actual computation. Reflective boundary conditions {dC/dz = 0) were imposed at the disk 
midplane due to the assumed symmetry. At the inner and outer radial boundaries, C=0 
was imposed throughout the simulations. At the inner boundary, this allowed materials that 
were migrating inwards to fall onto the star, being lost from the computational grid. This 
also meant that particles which were transported beyond the outermost radial grid were lost 
from the simulation. As the radial grid extended to ~80 AU this meant only a minimal 
amount of material was lost during a given simulation, and as the main area of interest here 
is for r < 30AU, this is expected to have minimal impact on the results presented here. 

A value of a;=10~'^ was assumed throughout the disk, with the viscosity given by v=acH, 
and the diffusivity of sohds equal to D = p/il+tlVt^) (Youdin and Lithwick 2007). In 



adopting this value of a, it is explicitly assumed that it is constant with time and location. 
However, if the MRI were responsible for driving disk evolution, it would operate most 
effectively in the regions of the disk in which ions were abundant enough for the disk to 
couple to magnetic fields. This would only be achieved in parts of the disk where significant 
ionization had occurred, meaning that parts of the disk would be located in the so-called 



"Dead Zone" (Gammie 1996 Glassgold et al. 1997). As the MRI would not operate here, it 



was thought such a region would not be subjected to the stresses that drive mass and angular 
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transport in protoplanetary disks, which would thus result in a different evolution than that 
described here. However, Fleming and Stone (2003) showed that overshoots of turbulent 
eddies from an MRI-active region could "stir up" stresses and turbulence in a Dead Zone, 



meaning that region would still be dynamically active. Further, recent work by Turner and 



Sano (2008) showed that viscous stresses could be produce in a neutral region of the disk. 



resulting in laminar (or less turbulent) accretion flows. This means flow structures similar 
to those studied in this work could still develop. The detailed impact of variations in a will 
be evaluated in future work. 



5. Model Results 



Cases of both outward and inward transport were examined. The outward transport 
of particles from the inner nebula was considered by looking at how materials formed at 
temperatures in excess of 1100 K were redistributed throughout the nebula. This temper- 
ature was chosen so that this simulation would be analogous to looking at how amorphous 
precursors could have been annealed in the hot inner disk and then transported outward. 



potentially to where comets formed in our own solar nebula (Nuth et al. 2000; Gail 2001 



Bockelee-Morvan et al. 2002 Ciesla 2007b), or to the large distances from the central star 
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Apai et al. 
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Wat- 



son et al. 2007). A higher temperature may be necessary if crystalline grains were formed 



purely from condensation (Wooden et al. 2007, and references therein). This would likely 
lead to slightly lower concentrations of materials in the outer disk as a smaller volume of the 
disk would be able to produce such high temperatures. However, given the steep thermal 
profiles in the disk (Figure 1), the difference in volume would not be large, so the change in 
concentrations would likely not be significant. 



For these outward transport calculations, a similar approach as Keller and Gail (2004) is 
adopted where the source of the materials is defined to be constant throughout the simulation 
(C = Co) where T > 1100 K and where z < H. This is similar to the approach used in 



Keller and Gail (2004), however, those authors considered the species as purely a vapor, so 



vertical settling or gas drag migration did not operate on their tracers. 

Inward transport of dust is also considered where materials that originate in the outer 
disk are tracked as they are redistributed throughout the protoplanetary disk. In this case, 
C = Co at the beginning of the simulation, but then the concentration evolves dynamically 
rather than being constant throughout the simulation. Such a situation would be akin to 
following the motions of short-lived radionuclides that are injected into the solar nebula from 



a nearby supernova (Ouellette et al. 2007), materials with anomalous oxygen isotope ratios 
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produced by CO self-shielding in the disk (Lyons and Young 2005), or material that rains 



onto the nebula after potentially being launched from an X-wind ( Shu et a7!]|1996 ). It should 
be noted that while this consideration is labeled "Inward Transport," outward transport of 
this material also occurs. It is so named because the general interest in such material or the 
natural application of these results has to do with how these materials would be mixed into 
the meteorite parent bodies and planets that form in the inner solar system. In the case of 
the "Inward Transport" calculations, the starting material originates in a radial span of ~10 
AU, centered at ~25 AU, and at a vertical span of one scale-height {H), centered at 2H at 
each radial location. 

In the case of the 10~^ Mq/jt disks, the models were run to a model time of 10^ years, 
during which 0.1 Mq of material would have been accreted onto the star. As this corresponds 
to roughly the mass expected for a typical protoplanetary disk around a solar mass star, this 
represented a logical stopping point. Similarly for the 10~^ Mq/jt disks, models were run 
to a model time of 10^ years for the same reason. The 10~^ Mq/jt were also modeled for 
10^ years (representing a period when 0.01 Mq was accreted by the central star) as typical 



protoplanetary disk lifetimes are on this order (Haisch et al. 2001). Also, in the case of our 



solar nebula, giant planet formation would likely have been at an advanced stage at that 



point (Hubickyj et al. 2005), meaning that the assumptions of axial symmetry and smooth 



surface density profiles likely would no longer hold for longer times. 



5.1. Outward Transport 

5.1.1. M=10-^ Mq/vt 

Figure 4 shows the temporal evolution of the concentration of solids formed at T >1100 
K for the case of M=10^^ Mq/jt and a=5 /im. Each contour in the plot represents a change 
in concentration by a factor of 2 compared to neighboring contours, with bright contours 
representing high concentrations and lower concentrations becoming progressively darker 
(the lowest contour has a value of 0.001 of Cq and the highest contour has a value of 0.512 of 
Co). With time, it is seen that outward transport does occur, with the concentration at the 



midplane exceeding the value at higher altitudes as reported by Ciesla (2007b). This is due 
to outward transport being favored along the midplane as this is the location of the outward, 
meridonal flows in the disk. In these flows, solids diffusing outwards are moving with the 
flow of the gas, rather than against it as would be the case in the accretional (towards the 
star) flows at higher altitudes or the net inward flows considered in ID models. 



Little difference is seen in these results and those where the outward transport of solids 
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with a=0.5 mm in a disk with M=10"® Mq/jt (Figure 5), despite these latter particles 
being two orders of magnitude larger in size (and six orders of magnitude larger in mass). 
This is due to the high temperatures and high gas densities that are present in a young 
protoplanetary disk which can accrete at a rate of M=10~® Mq/jt. Both of these factors 
contribute to reducing the stopping time of the dust particles and making them very small 
compared to the orbital timescales of the nebula at all heights. As a result, the two particle 
sizes have very similar dynamic responses due to their interactions with the gas, meaning 
the distributions of solids at the end of the simulations are nearly the same. Thus, in both 
cases particles are carried well beyond 20 AU after 10^ years of evolution. 



5.1.2. M=10-^ Mq/vt 

Particles are not carried outwards as rapidly in the case of M=10~'^ Mq/jt (Figures 6 
and 7) as the outward flows are lower in magnitude at the midplane and because the 1100 
K isotherm is located closer to the star. However, even after just 10^ years of evolution, 
materials from the hottest regions of the nebula still are transported beyond 20 AU, result- 
ing in concentrations of a few percent in the region where comets like Wild 2 would have 



originated (Charnoz and Morbidelli 2007). Here, though, because the disk is less massive 
and cooler, contrasts in the dynamics of the different sized particles begin to be more readily 
seen. Firstly, the larger particles are not lofted as high as the smaller dust particles as shown 
at the 100,000 year mark, where outside of 15 AU small dust particles are present 4 AU from 
the disk midplane, whereas the greater settling velocities of a=0.5 mm particles keep them 
more concentrated around the midplane. As the larger particles are kept at lower altitudes 
as compared to the smaller dust, they spend longer periods of time in the regions where 
outward transport occurs more efficiently. As a result, the concentration of the a=0.5 mm 
particles around the midplane in the outer disk slightly exceeds that of the a=5 /xm particles. 
This effect is most noticeable in this region because the lower gas densities and the lower 
temperatures produce more significant differences in the stopping times of the particles than 
would be found in the hotter, denser, inner nebula. 



5.1.3. M=10-^ Mq/vt 

Figures 8 and 9 show the results for the M=10~^ Mq/jt cases. Again, the drop in 
accretion rate results in lower outward flux of high temperature materials due to the lower 
velocities at the disk midplane and due to the fact that high temperatures are limited to the 
very inner edge of the nebula. This decrease in outward flux means that particles of either 
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size make it to just beyond 10 AU in 100,000 years, rather than beyond 20 AU as in the 
higher mass accretion rate cases. Again the lower temperatures in the disk and lower gas 
densities result in greater differences in the dynamics of the two sizes of particles. In this 
case, even in the inner disk, the a=0.5 mm particles are not lofted to the same heights as the 
a=5 /xm dust. As before, this also results in higher concentrations of the larger particles in 
the outer parts of the disk (at 10 AU in the 100,000 year mark, for example). This is a result 
that is counter to what would be found in the ID models typically used in the past, as the 
greater inward velocities of the larger particles due to gas drag and lower diffusivities of the 
particles would hinder their outward transport as compared to their smaller counterparts. 



5.2. Inward Transport 

5.2.1. M=10-^ Mq/vt 

Here the "Inward Transport" cases are described, where the dynamics of particles origi- 
nating at radial distances of 20-30 AU and 1.5-2.5 scale heights above the disk midplane are 
calculated. Figure 10 shows the results of the case of a=5 /xm and M=10^^ Mq/yt. The 
contours in these figures are the same as those used in the previous section. As expected, 
these materials are redistributed throughout the disk, with inward transport occurring most 
rapidly along the disk surface, where the inward flows of gas due to viscous evolution are 
the greatest. Thus immediately after these grains are "injected" or "created" (depending on 
what kind of materials they represent) in the disk, there are strong vertical gradients in the 
concentration of these species. These gradients persist throughout the simulation, even as 
diffusion works to smooth them out. This is due to the fact that while dust mixes down to 
the midplane due to gravity and diffusion, this process lags behind the rapid inward motions 
of the grains at higher altitudes. 

Gravitational settling occurs more rapidly for the 0.5 mm particles as shown in Figure 
11. While vertical gradients similar to the type seen in the a=5 //m case are present during 
the early part of the simulation, after ~50,000 years, the 0.5 mm particles cluster around the 
midplane where their vertical distribution approaches the canonical steady state situation- 
with vertical diffusion balancing vertical settling. (The steady state is approached but not 
fully reached because the vertical adjustment is not instantaneous, and in the presence 
of radial transport, the conditions throughout a vertical column are constantly changing.) 
Thus, vertical gradients are also present, but with the highest concentrations of 0.5 mm 
particles being achieved at the disk midplane and decreasing with height, as opposed to the 
other way around for the smaller particles. 
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That the a— 5 /xm grains remain lofted at high altitudes, and thus are entrained in the 
rapid inward flows of the disk, for longer periods of time make them easier to deliver to 
the inner disk. This is shown in detail in the different panels in Figure 12. Panels A and 
B show the concentrations of the 5 fim. and 0.5 mm solids in the inner disk at the end of 
the simulation respectively. It can be seen that the contours extend to smaller heliocentric 
distances in the case of the smaller dust particles than the larger dust particles. This is 
reinforced in the plot in panel C, which shows the concentrations of both sized particles 
along the disk midplane at the end of the simulation. Inside of 10 AU, the concentration 
of the smaller particles exceeds that of the larger ones everywhere despite the more rapid 
settling velocities of the larger particles. This also means that the small particles are accreted 
onto the star, and thus lost from the disk, more readily than the larger particles. This is 
shown in panel D, where the fraction of the injected material remaining in the disk is plotted 
as a function of time. The fraction of small dust decreases more rapidly than the larger dust 
particles. While the fraction remaining in the disk at the end of the simulation is still large 
(~95%), this reflects the fact that the materials of interest started so far from the inner edge 
of the disk (>20 AU). Had this material instead been injected at 10 AU, the fraction lost, 
and discrepancy between the two different sized particles would likely have been larger. 



5.2.2. M^lOr^ M^/yr 

Figure 13 shows how the small dust grains are transported over time in the M=10~^ 
Mq/yy disk. The overall evolution of the solids again matches general expectations in that 
the solids begin to settle out of the surface layers, while being redistributed radially due 
to the large-scale flows associated with viscous evolution and the diffusion caused by the 
turbulence. While the particles do get redistributed throughout the disk, with outward 
transport carrying them all the way to the outer boundary of the disk, their motions are 
largely inward, as is expected given the net movement of disk materials toward the star under 
viscous evolution. Thus after 1 million years of evolution, the particles are concentrated to 
the highest levels inside of 10 AU. 

Note again that vertical gradients persist immediately after injection onto the surface 
of the outer disk, though their magnitude at 100,000 years is not as great as in the case of 
M=10~® Mq/jt. The vertical diffusion timescale for these particles, say to make it down 
to the midplane from where they were injected, is (2iJ)^/z/=4/7^/«ciJ=4(aQ/^)^^, meaning 
it is independent of the disk structure or mass accretion rate. The radial flows in a disk, 
however, decrease in magnitude as the mass accretion rates drop, meaning that the conditions 
in a column do not vary as rapidly, allowing vertical diffusion to smooth out gradients more 
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effectively. 

Figure 14 shows the results for the M—10~'^ Mq disk, but with the larger particles. 
Again, the large-scale flows and turbulence redistribute the particles throughout the disk. 
However, because of the lower densities and temperatures in this disk as opposed to the 
M—10~^ Mq/jt case, the larger particles settle towards the midplane much more rapidly, 
resulting in very different dynamical evolution when compared to the smaller particles. The 
more rapid settling allows larger particles to quickly cluster in the region of the nebula where 
outward transport is favored (around the midplane), slowing their delivery to the inner disk 
and resulting in greater mass fluxes to the outer disk. This is more clearly seen in Figure 
15, in which the concentrations of the two different species in the inner disk are shown in 
the top panels at the end of the simulation. There are not significant differences in terms of 
the concentrations in the inner part of the disk (panels A and B), however, panel C shows 
that the concentration of the larger species at the midplane well beyond 40 AU exceeds the 
concentration of the smaller dust particles. This is in part due to the greater settling rates of 
the larger particles, but also due to the fact that the outward flows in this region preserve the 
larger particles from migrating to the inner edge of the disk as rapidly as the small particles. 
This is reinforced in panel D, which again shows that the smaller particles are lost from the 
disk more rapidly than the larger particles, with the differences being larger here than in the 
M=10-^ Mq/yt case. 

5.2.3. M^lOr^ M^/yr 

Figures 16 through 18 show the results for a simulation with M=10~^ Mq/jt. Again, 
this disk has lower densities and temperatures than those disks considered above. As a 
result, the large-scale ffows have decreased in magnitude, leading to lower levels of radial 
transport. Vertical gradients in the concentration of materials are minor, with the small 
particles becoming well mixed in the vertical direction and the large particles clustering 
around the midplane. The outward flows around the midplane decrease but are still present; 
however the lower velocities mean that they are not as effective at slowing the inward move- 
ment of the larger particles due to diffusion or gas drag. As a result, there is little difference 
between the delivery of the materials to the inner disk. Panel C of Figure 18 shows that the 
concentration of the large particles at the disk midplane exceeds that of the small particles 
inside of ~40 AU. This is largely due to the greater settling efficiency of the larger particles 
and not due to more rapid inward transport of the small particles in this case. That the 
smaller particles have greater abundances at the midplane beyond 40 AU is due to the fact 
that gas drag migration serves to prevent the larger particles from migrating outwards to 
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the same levels as in the previous cases. And while the greater settling efficiency of the large 
particles removes them from the rapid inward flows that occur along the disk surface, the 
difference in the fraction of material lost to the star due to the accretion process is minor 
for the different sized particles (Panel D), though again the larger particles are still retained 
more efficiently than the small particles. 



6. Summary and Discussion 



The results of the different models presented here demonstrate that the vertical varia- 
tions in the rates and preferred directions of radial transport of solids in a protoplanetary 
disk are important to consider when talking about the early evolution of primitive materials. 
These variations arise due to the fact that the gas densities, pressure, and their respective 
gradients, all of which are responsible for driving motions in a protoplanetary disk, are strong 
functions of location. In terms of the viscous evolution of the disk, the vertical variations in 
gas density are responsible for driving a flow structure in which the gas moves toward the 
star rapidly at high altitudes in the disk, but moves outward, away from the star around the 
disk midplane. The inward motions of the gas are responsible for the mass accretion onto the 
central star, representing its final stages of pre-main sequence growth, whereas the outward 
motions are responsible for the outward transport of angular momentum to compensate. 
The relative fraction of material transported outward at a given location can range from ~5 



to 30% of the mass that is transported inward at a given time (Keller and Gail 2004). 



Further, the gas drag induced velocities of solids in the solar nebula will also vary 
with height as the radial pressure gradient, which is responsible for the gas rotating at a 
non-Keplerian rate, also varies with height. Unlike the large-scale flows caused by viscous 
evolution, however, the gas drag induced velocities will reach their largest, negative (towards 
the central star) velocities near the disk midplane. At higher altitudes, the inward velocities 
diminish, and in fact, will become positive as the radial pressure gradient switches signs, 
leading the gas to orbit at super-Keplerian rates at higher altitudes. 

These variations in the radial transport models are overlooked in traditional one-dimensional 
models of material transport in protoplanetary disks (e.g. Stepinski and Valageas|[l997 Gail 



2001 Bockelee-Morvan et al. 2002; Cuzzi et al. 2003 Ciesla and Cuzzi 2006). Instead, the 



motions of the particles due to the viscous evolution of the disk in these models are found 
by assuming that those particles are uniformly suspended throughout the height of the disk 
(thus ignoring vertical settling) and their net velocity is found by calculating the vertically 
averaged motions of the gas at that location (l^=-3z//2r under steady-state condtions). For 



those ID models that account for the radial motions due to gas drag (e.g. Stepinski and 
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Valaereas |[l997l [Cuzzi et a/.|[2003| [Ciesla and C'uiil|[2006| , these velocit ies are generally found 



by calculating the inward velocity at the disk midplane and applying that velocity to all 
particles, regardless of where they are located in the vertical column above the disk mid- 
plane. While this may be appropriate for particles which are tightly clustered around the 
disk midplane due to vertical settling (~meter-sized bodies which do not experience signifi- 
cant vertical diffusion, for example), this would overestimate the inward gas drag velocities 
of particles lofted to higher altitudes. Thus while ID models are instructive for investigating 
the general dynamical evolution of materials in protoplanetary disks, consideration of the 2D 
effects can lead to important differences in the dynamical evolution of primitive materials. 



Here, as in Ciesla (2007b), it was shown that the outward flows that are generated 
around the midplane in a viscously evolving disk aid in the outward transport of materials 
by providing an environment in which materials can diffuse outwards without being hindered 
by the accretional flows associated with protoplanetary disk evolution. This provides a path 
by which high temperature materials (refractory inclusions and crystalline grains) collected 
from comet Wild 2 or seen in comet Halley and Tempel 1 could be delivered to the outer 
solar nebula. Predicting the particular fraction of silicate grains in the outer disk that would 
be crystalline in each of the models presented here is not straightforward, unfortunately. 
This fraction depends on the amount of material that is delivered to the outer disk as well 
as the amount of unprocessed material that is present. The abundance of this material will 
be determined by its dynamical interaction with the gas in the outer disk, the extent of the 
disk (the initial inventory of unprocessed materials) and the rate at which these unprocessed 
materials coagulate into larger bodies. A first order approximation, however, would be that 
the total processed and unprocessed silicate grain abundance would add up to the solar ratio 
{C crystal + C amorphous ~ 0.005). This wouM suggcst that the crystalline grain fraction at 10 
AU at the end of each simulation would be ~40% for each of the M = lO~^M0/yr cases, 
~10% for the M = lO"'^M0/yr cases, and < 1% for the M = lO"^M0/yr cases. These 
numbers should only be taken as rough estimates, however, and future work should follow 
both the unprocessed and processed materials simultaneously. 



These numbers are generally comparable to those found by Wehrstedt and Gail (2008) 



in their two-dimensional models, though direct comparisons are not straightforward due to 



different model treatments. For example, Wehrstedt and Gail (2008) assumed that grains 



could become crystalline at temperatures as low as 800 K, allowing for a larger volume of 
grains to be processed at any given time. More importantly, however, as discussed above, 
these authors focused on the dynamics of dust particles in the region around the midplane, 
rather than considering the dynamics of solid particles at very high altitudes. These higher 
altitudes still contribute significantly to the inward mass fiux in the disk. Thus these models 
likely overestimate the outward fiux of high temperature materials and underestimate the 
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inward flux of unprocessed materials. It should be noted that the differences due to the 



particular choice of viscosity model (a-viscosity in this work, /5-viscosity in Wehrstedt and 



Gail (2008)) needs to also be quantified in future studies. 



Not only does the rate of mass accretion in the disk determine the amount of outward 
transport that occurs in a disk, but the timing of this transport as well. For disks with mass 
accretion rates in excess of M >1Q~'^ Mq/ji, high-temperature materials can be transported 
beyond 20 AU on timescales of ~10^ years. Such rapid outward transport is likely necessary 
given that the formation timescale of the giant planets via core accretion is approximately 
1-2 Myr ( Hubickyj et a/.||2005 ). Delivery of high temperature materials on longer timescales. 



such as a few million years, would be problematic as the outward transport of the grains 
would likely be blocked by the young giant planets, either because the cores would accrete 
the grains themselves or because the planets would begin to open a gap in the disk through 
which particle transfer would likely be difficult. 

At lower mass accretion rates, with M ~lO~®M0/yr, typical of 10^ year old T-Tauri 



stars (Hartmann et al. 1998), delivery of crystalline silicates or CAIs to the outer disk is 
difficult as high temperatures are generated only at the very inner edge of the disk. The 
lower outward flow velocities along the disk midplane and the smaller volume of the disk 
in which high temperature materials can be formed or processed means that delivery of 
crystalline grains and refractory inclusions to beyond 20 AU will not occur as quickly as in 
the higher mass accretion rate models. Thus it is likely that the high temperature material 
seen in comets were processed and delivered to the outer solar nebula during the early stages 
of disk evolution. 

However, even at lower mass accretion rates, such as M ~lO~^M0/yr, the flows at the 
disk midplane are likely still important in terms of understanding the dynamical evolution 
of sub-cm sized objects, as the outward flows velocities in the inner disk are >10 cm/s in the 
model presented here. For the a=0.5 mm particles (which is roughly the size of chondrules) 
considered here, this means that the outward flows at the disk midplane exceed the inward 
drift velocities at the disk midplane that arise due to gas drag (Figure 3) out to about 20 
AU. In fact, given that inward drift velocity for small particles is roughly proportional to 



its stopping time ( Cuzzi and Weidenschilling 2006 ) the outward flows at the disk midplane 



would exceed the inward drift velocities of objects with radii of 5 mm (comparable to the 
larger CAIs found in CV chondrites) out to a few AU, offering a way of preserving these 
objects in the nebula for millions of years before the formation of chondrules and accretion 



of chondrites (Amelin et al. 2002 Kita et al. 2005). 



The variations in flow structure also are important to consider when discussing the 
inward transport of materials from the outer parts of the protoplanetary disk. Materials 
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that would be injected into the solar nebula, rain down after being lofted upward in a 
protostellar jet, or are created due to photochemistry on the disk surface would initially be 
subjected to rapid inward flows in the very upper regions of the disk. This would remain 
true until the particles are transported down towards the midplane where the inward flows 
diminish in magnitude or the flows reverse direction. This removal from the upper part 
of the disk occurs more rapidly for larger particles due to their greater settling velocities, 
meaning that smaller particles can be transported to the inner disk on shorter timescales 
than larger ones. In fact, because smaller particles can remain lofted at high altitudes for 
long periods of time, they are more likely to be accreted onto the central star in accretionary 
flows as compared to the larger bodies that settle around the disk midplane as demonstrated 
in the simulations presented above. This remains true as long as the velocities due to viscous 
evolution of the disk exceed those of that arise due to gas drag. This result is not limited to 
materials that begin at high altitudes. Vertical diffusion would loft small particles to higher 
altitudes than larger ones, thus introducing them to the regions of rapid inward flow more 
readily. This runs counter to the previous ideas that millimeter to centimeter sized objects 
would be more difficult than fine dust to retain in the solar nebula due to their inward drift 
velocities caused by gas drag. Instead, it appears that aggretation into bodies of this size 
may actually have allowed them to be retained in the solar nebula more effectively. 

It should be noted that millimeter-sized grains themselves do not necessarily serve as 
good analogs for the materials that would be injected into the solar nebula by a supernova 
dOuellette et a/.|2007D or for the heavy {^^O -'3oor) water ice that would be created by CO self- 
shielding in the outer nebula ( Lyons and Young||2005 ) , as these materials are likely contained 
within or adsorbed by grains much smaller and closer to the a=5 fim cases discussed above. 
However, it should be kept in mind that these smaller grains would likely have experienced 



collisions with other grains that led to sticking and growth. Weidenschilling (1997) showed 
that micron-sized grains would coagulate into millimeter-sized grains at 30 AU on timescales 
of less than 10,000 years. Those coagulation calculations assumed a non-turbulent and non- 
evolving disk; growth at these small sizes may actually be accelerated by turbulence as the 



random motions of the dust particles would lead to more frequent collisions ( 


Weidenschilling 


1984 


DuUemond and Dominik 


2005 


Ciesla 


2007a 


Ormel and Cuzzi 


2007 


). Growth to 



very large sizes may be inhibited due to turbulence, as larger particles will develop relative 
velocities that become destructive in collisions, rather than accretional, but millimeter-sized 



grains appear to be below this threshold (see discussion in Cuzzi and Weidenschilling 2006). 
Thus, the results of the a=0.5 mm models are likely more reflective of how solids from the 
upper layers of the outer disk get transported to the inner disk. 

The issue of retention is particularly important for the case of injection of short-lived 



radionuclides directly into the solar nebula (Ouellette et al. 2007). If this material was in- 
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jected, it would likely have been done in one or more Rayleigh- Taylor fingers (Boss 2007), 



with materials being peppered along the surface of the nebula as assumed here. Ouellette 



et al. (2007) demonstrated the ratios of ^^Al/^^Al and ^°Fe/^^Fe would be close to the in- 
ferred meteoritic values if the solar nebula were located 0.15 pc from a supernova. These 
calculations, however, assumed that these materials became well mixed in the nebula. How- 
ever, as shown here, injected material that was present as small grains (<10 /zm, as expected 
for the supernova grains) at the surface of the nebula would be lost to the Sun faster than 
larger solids, which would have settled around the midplane. Thus to ensure that these 
injected materials become well mixed in the solids, they would quickly have to be accreted 
into the larger bodies that had already formed. If they were not, the injected material would 
be lost at a faster rate than the native solids, thus leading to lower ratios of the abundance 
of short-lived radionuclides to stable isotopes. This would be particularly true if solids were 
injected closer to the Sun than the ~20 AU distance considered here as accretion onto the 
sun would happen on shorter timescales. This would require the solar nebula to be closer 
to the supernova to achieve the meteoritic ratios. Such considerations should be made in 
future studies. 

An important point to note in each of these simulations is that the concentration at the 
end of each of the "Inward Transport" simulations shown here is not uniform throughout 
the disk, but rather, gradients exist. This is particularly seen in panel C of Figures 12, 15, 
and 18, which show that the concentrations of the species at the disk midplanes vary with 
distance from the star. Radial mixing timescales for cases such as these are typically defined 
as the diffusive timescale, tmix=i^'^ / ^ ■, which for the inward transport calculations presented 
here would be 6x10^, 1.5x10^, and 2.5x10^ years at 20 AU for each of the disks considered. 
These numbers represent estimates for the amount of time it would take for the "injected" 
material tracked here to mix into the very inner part of the protoplanetary disk. While the 
model simulations are not run to these times for the cases considered here, it can still be 
seen that delivery of significant amounts of material occurred on shorter timescales. However, 
complete homogenization (little to no concentration gradients) of the injected material had 
not been achieved at the end of the simulations, nor is it expected to be achieved at any 
reasonable time afterward. The evidence for this is that gradients exist at every point in 
the disk, regardless of the distance from where the materials originated. Because of the 
variation in the direction and magnitudes of the large-scale flows in the disk, conditions 
are constantly changing, not allowing diffusion to work to smooth out the concentration 
gradients completely. 



The agreement between Pb-Pb ages and Al-Mg ages of CAIs and chondrules (Amelin 



et al. 2002) suggests that ^^Al was uniformly distributed relative to the stable isotope ^^Al 



throughout the inner solar nebula where these materials formed, implying that if injected. 
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^^Al was homogeneously distributed throughout the inner solar system where these objects 
formed. Such apparent homogenization in the injection model could be achieved if mixing 
had occurred in such a way that these objects formed near the concentration maximum, 
where the radial gradients in concentration are the lowest (but non-zero). Alternatively, 
near homogenization of injected materials could have been achieved if the materials were 
contained in small dust grains and they were injected at a late stage of nebular evolution. 
Radial concentration gradients are smallest in the cases presented for the 5 /xm grains in the 
disk accreting at lO^^Mo/year, as the large-scale flows slowed to the point that diffusion 
became more effective at smoothing out gradients. However, this low mass accretion rate 
also would not produce very high temperatures in the disk, meaning that CAIs are unlikely 
to form or their formation and redistribution was limited to only regions very close to the 
Sun. Further, the small dust grains would likely have coagulated into larger objects, such as 
the 0.5 mm grains considered here, which show comparable or stronger radial gradients in 
all cases considered. 

Therefore, if ^^Al was injected directly into the nebula, it would likely not have been 
perfectly homogenized, but rather gradients, or variations in its concentration, would have 
existed. These variations would lead to apparent age differences between chondritic com- 
ponents, the magnitude of which depended on the particular values of ^^Al/^^Al that were 
recorded by those components. For example, a 10% variation in the ^^Al/^^Al ratio could 
easily develop over a scale of a few AU based on the results of the simulations presented 



here, and also is in line with the variations found by Boss (2008) for materials injected 



into a marginally gravitationally unstable disk. If two solids formed contemporaneously, 
one with ^^A1/^''A1 =5x10^^, the canonical value within CAIs (see Jacobsen et al. 2008, for 



the most recent discussion on this issue), and one with ^^Al/^^Al =4.5x10"^, this would be 
interpreted as an age difference of ~100,000 years, which is less than the uncertainty on the 
Pb-Pb ages. If this were the case, and ^^Al was injected directly into the solar nebula, the 
inherent uncertainty in Al-Mg ages would be at least ~100,000 years. 



Recently, measurements of Fe, whose presence in solar system materials (e.g. Tachibana 



and Huss 2003) has been among the motivating factors for looking at the direct injection 



of short-lived nuclides into the young solar system from a supernova, have helped to shed 
light on how well homogenized potential supernova derived materials were in the solar neb- 
ula. Bizzarro et al. (2007) originally reported measurements of ^°Ni deficits in meteorites 
and terrestrial samples that indicated ^°Fe was not present in planetesimals that accreted 
within the first 10^ years of the formation of the solar system. This finding supported the 
model of direct injection into the solar nebula by a supernova dOuellette erar]|2007D . In this 
scenario, the 10^ year interval represented the time it took for ^*^Fe to be injected and mixed 



in to the inner solar nebula. More recently however, Dauphas et al. (2008) analyzed similar 
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samples as Bizzarro et al. (|2007l), but found that ^°Fe was homogenized in the terrestrial 



planet region of the solar nebula such that its concentration varied by less than 10%. Again, 
these results could be achieved in the injection model if the maximum in the concentration 
of ^°Fe had migrated to the inner solar system at the time that planetesimals there began to 
form. If this were not the case, or if the variation in ^°Fe is much less than 10% in planetary 
materials, the origin of short-lived radionuclides may be best understood as resulting from 



the injection into the natal molecular cloud prior to (or as the cause of) collapse (Boss et al. 



2008). 



The results presented here differ somewhat from those of Boss (2008) who showed that 



the concentration of "dye" that was injected to a marginally gravitationally unstable disk 
would be mixed to the point that there would be less than 10% variation in its concentration 
throughout the entire disk after ~1000 years. This may represent a fundamental difference in 
how transport occurs in gravitationally unstable disks and viscous disks, or it may represent 
differences in the ways that calculations were performed. For example, while the study 



described here tracked the evolution of solids, the dye in Boss (2008) was more akin to vapor 



in the disk, meaning that the dynamics of vertical settling and gas drag were not considered. 
While this may approximate the dynamical evolution of small dust particles which are well 



coupled to the gas, as described by Haghighipour and Boss (2003), these processes are likely 



important as solids grow beyond a few microns in size. Future work aimed at determining 
the differences in how solids evolve in marginally gravitationally unstable and viscous disks 
will help to identify which model better describes the early evolution of our solar nebula. 
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Model Parameters for Star-Disk System 



Stellar Mass 




Stellar Temperature 


4000 K 


Stellar Radii 


3 Rq 


Grazing Angle, 


0.05 


Disk Opacity, k. 


5 cm^/g 


Disk Inner Radius, rj„ 


0.05 AU 


Disk Outer Radius 


84 AU 



Structural Properties of Disks Considered 



M (Mo/yr) 


10-6 


10-7 


10-8 


r-iioo K (AU) 


3.4 


1.3 


0.5 


at 1 AU (AU) 


0.14 


0.088 


0.056 


H at 20 AU (AU) 


3.0 


2.0 


1.5 


Time Computed (yrs) 


10^ 


10^ 


10^ 
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Fig. 1. — The surface density and temperature profiles for the protoplanetary disks consid- 
ered here with mass accretion rates of 71^=10^*^, 10^^, and lO^^M0/year. 



- 36 - 




Fig. 2. — The velocity profile for the gas in an viscously evolving protoplanetary disk. 
Arrows indicate the direction of the flow, with the length of the arrow being proportional to 
the magnitude of the flow. For reference, the outward flow rates of the gas at the midplane 
for the different disks considered here are given in Figure 3. 
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Fig. 3. — The outward flow velocities of the gas (left) and the inward velocities of the 0.5 
mm particles due to gas drag (right) both along the disk midplane. Gas drag velocities for 
the 5 fim grains would be ~100x smaller. Outward flow velocities decrease with decreasing 
mass accretion rates, while gas drag velocities increase in magnitude. As can be seen, only 
in the outer regions of the M=lO^^M0/year disk is the net velocity of the 0.5 mm particles 
inward at the disk midplane. 
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Fig. 4. — Time evolution of the concentration of high temperature 5 /im grains in the pro- 
toplanetary disk with M=lO~^M0/year. The concentration of materials was held constant 
at Co=l wherever T>1100 K and z < H. Contours represent changes in concentration by a 
factor of 2, ranging from 0.001 to 0.512. 
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Fig. 5. — Same as Figure 4, with a 



,=0.5 mm. 
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Fig. 7. — Same as Figure 6, with a— 0.5 mm. 
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Fig. 8 — Same as Figure 4, with M^lO^^Mo/year. 
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Fig. 9. — Same as Figure 8, with a— 0.5 mm. 
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Fig. 10. — Time evolution of the concentration of 5 //m grains in the protoplanetary disk 
with M =10^^ Mq / jeai . The concentration of materials was set at Co=l at t=0 between 20 
and 30 AU and at heights between 1.5 to 2.5 H. Contours are the same as those in Figure 
4. 
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Fig. 11.- 



Same as Figure 10, with a=0.5 mm. 
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Fig. 12. — A summary of the results of the inward transport models for the disk with 
M=lO^®M0/year. Panels A and B show the distribution of the 5 /^.m and 0.5 mm grains 
at the end of the simulation respectively, with the contours having the same meaning as in 
Figure 4. Panel C shows the concentration of the different sized particles particles at the disk 
midplane at the end of the simulation. Radial gradients in the concentrations are apparent. 
Panel D shows the mass fraction of solids retained in the disk as a function of time. That 
material that stays lofted at high altitudes and is pushed inward by the viscous flows of 
the disk begins to be accreted onto the central star after ~10^ years. Because larger grains 
settle out of these upper regions more readily, they are retained more effectively than the 
smaller grains, an effect that runs counter to expectations based on previous one-dimensional 
models. 
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Fig. 13. — Time evolution for the 5 fim grains for the M—10 ^M0/year disk. The initial 
concentration was defined as in Figure 10. 





Fig. 14. — Same as Figure 12, with a— 0.5 mm. 
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Fig. 15. — Same as Figure 11, with M=10 ''M0/year. Again, radial gradients remain in the 
disk everywhere and the smaUer grains are lost from the disk faster than the large grains. 
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Fig. 17, 



, — Same as Figure 14, with a— 0.5 mm 



-49- 





Fig. 18. — Same as Figure 11, with M=10~^MQ/year. The gradients in the distribution of 
the 5 iim grains are the smallest of the three simulations, as the large-scale flow velocities 
had diminished so that they were not significantly altering the gradients that diffusion was 
working to smooth. This was not the case for the 0.5 mm grains, as gas drag velocities were 
high enough to constantly change the conditions on which diffusion was operating. 



